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ABSTRACT 


Theoretical expressions are derived for the steady-state frequency response of 
a thermocouple wire. In particular, the effects of axial heat conduction are 
demonstrated for a non-uniform wire with unequal material properties and wire 
diameters across the junction. The amplitude ratio at low frequency co -» 0 agrees 
with the results of Scadron and Warshawsky (1952) for a steady-state temperature 
distribution. Moreover, the frequency response for a non-uniform wire in the 
limit of infinite length / _ > oo is shown to reduce to a simple expression that is 
analogous to the classic first order solution for a thermocouple wire with uniform 
properties. 

Theoretical expressions are also derived for the steady-state frequency 
response of a supported thermocouple wire. In particular, the effects of axial heat 
conduction are demonstrated for both a supported one material wire and a two 
material wire with unequal material properties across the junction. For the case of 
a one material supported wire, an exact solution is derived which compares 
favorably with an approximate expression that only matches temperatures at the 
support junction. Moreover, for the case of a two material supported wire, an 
analytical expression is derived that closely correlates numerical results. 

Experimental measurements are made for the steady-state frequency response 
of a supported thermocouple wire. In particular, the effects of axial heat 
conduction are demonstrated for both a supported one material wire (type K) and a 
two material wire (type T) with unequal material properties across the junction. 
The data for the amplitude ratio and phase angle are correlated to within 10% with 
the theoretical predictions of Forney and Fralick (1991). This is accomplished by 
choosing a natural frequency co n for the wire data to correlate the first order 
response at large gas temperature frequencies. It is found that a large bead size, 
however, will increase the amplitude ratio at low frequencies but decrease the 
natural frequency of the wire. The phase angle data are also distorted for imperfect 
junctions. 
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1.0 UNSUPPORTED THERMOCOUPLE WIRE: THEORY 


Theoretical expressions are derived for the steady-state frequency response 
of a thermocouple wire. In particular, the effects of axial heat conduction are 
demonstrated for a non-uniform wire with unequal material properties and 
wire diameters across the junction. The amplitude ratio at low frequency co — > o 
agrees with the results of Scadron and Warshawsky for a steady-state 
temperature distribution. Moreover, the frequency response for a non-uniform 
wire in the limit of infinite length /-><» is shown to reduce to a simple 
expression that is analogous to the classic first order solution for a thermocouple 
wire with uniform properties. 

1.1 INTRODUCTION 

The measurement of fluctuating temperatures in the high speed exhaust of 
a gas turbine engine combustor is required to characterize the local gas density 
gradients or convective heat transfer (Fralick, 1985). Although thermocouples 
are suitable for the measurement of high frequency temperature fluctuations (< 
1 KHz) in a flowing gas or liquid, the measured signal must be compensated 
since the frequency of the time dependent fluid temperature is normally much 
higher than the corner frequency of the thermocouple probe (Scadron and 
Warshawsky, 1952). Moreover, the amplitude and phase angle of the 
thermocouple response may be attenuated by axial heat conduction for rugged 
thermocouples of finite length (Elmore, et al.; 1983, 1986). 
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In the present chapter, the theoretical steady-state frequency response of a 
thermocouple wire has been calculated to include the effects of axial heat 
conduction. Solutions are derived for a non-uniform wire with unequal 
properties across the junction (e.g., different wire diameter, thermoconductivity, 
etc.). Solutions are presented in the form of the amplitude ratio and phase angle 
for both the uniform the non-uniform thermocouple wires. In both 
configurations, the results are presented for thermocouple wires of either finite 
length where the effects of axial heat conduction are demonstrated or wires of 
infinite length for which the neglect of axial conduction is justified. For the case 
of the uniform wire, the theroetical results are compared with the series solution 
of Scadron and Warshawsky (1952). Moreover, for non-uniform wires the 
results are compared in the limit of low frequency with the steady-state 
temperature solution of Scadron and Warshawsky. 


1.2 THEORY 

The steady-state frequency response of a thermocouple wire will be 
developed in the following sections with the following assumptions: (a) the 

amplitude of the fluctuating fluid temperature is small relative to the mean 
absolute temperature (b) the thermocouple dimensions are small relative to the 
size of the turbulent eddies or enclosure dimensions (c) radial temperature 
gradients in a wire cross section can be neglected and (d) radiative heat transfer 
can be neglected relative to conduction and convection. 
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1.2.1 Uniform Wire 


In this section the geometry of Fig. 1 is considered where the material 
properties of thermal conductivity k, specific heat c and wire density pare 
assumed to be equal on both sides of the thermocouple junction. If the probe is 
subjected to a flowing fluid, whose temperature fluctuates about its mean, the 
expression for the local conservation of energy in the thermocouple wire 
becomes (Scadron and Warshawsky, 1952) 


5T, a 2 T w 4h 


iw. = (T -T ) 

^ ~ ,2 ^ V A w' 


at " dx 2 pcD 8 -w ' d.i) 

where a = k/pc is the thermal diffusivity of the wire, T g is the ambient fluid 
temperature, h is the convective heat transfer coefficient, D is the wire diameter 
and T w is the local wire temperature measured along the axis at a distance x from 
the centerline of Fig. 1 . 

If the wire and fluid temperatures are measured relative to the support and 
mean fluid temperature To, the variation in the ambient fluid temperature can 
be written in the complex form 


T g (t) = T 0 +T f e 


icot 


( 1 . 2 ) 


where cj is the angular frequency of the ambient fluid. Since Eq. (1.1) is linear, 
we now seek a solution for the local wire temperature of the form (Hildebrand, 
1976) 


T W = T 0 + T m (x)e ia> '. (1-3) 

Referencing all temperatures with respect to the steady-state temperature T 0 and 
normalizing with respect to the amplitude of the fluctuating ambient fluid 
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temperature Tf, one defines a local normalized steady-state frequency response 
T(x ) for the thermocouple wire of the form 


t w -r, _ Tjxy 


= T(x)e im . 


(1.4) 


V v 

Substituting Eqs. (1.2) and (1.3) into Eq. (1), one obtains an ordinary 
differential equation of the form 


d*T 

dx 2 


icoT = a ^ 7T + co n (l-T) 


(1.5) 


where T = f(x) is the frequency response and the vector notation will be dropped 
for simplicity. Thus for the geometry of fig. 1, one seeks a solution to the non- 
homogeneous linear second order differential equation for the dependent 
variable T of the form 

y T" - G(co)T = - 1 (1.6) 

subject to the boundary conditions 

T(l) = T(— /) = 0. (1.7) 

Here, the parameters in Eqs. (1.5) and (1.6) are defined as 


co. = 


Ah 


a 

(o. 


G((o) = 1 + /(— ) 

co. 


(1.8) 


" pcD' 

where co n is the natural frequency of the uniform wire of fig. 1. 

The general solution to Eq. (1.6) can be written in the form (Hildebrand, 
1976) 


T = A sinh qx + B cosh qx + 


G((o) 


(1.9) 
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where the constants A and B are complex, 1 /G(q)) represents the particular 
solution and the parameter 


q = 



Substituting x - + 1 in Eq. (1.9), the boundary condition Eq. (1.7) yields values for 
the constants 


A = 0, B = 

G(co)coshql (1.10) 

Hence, one obtains a steady-state frequency response for the simple uniform wire 
of Fig. 1 in the form (Forney, 1988) 


T = 


1 ^ cosh qx 

G(co) cosh ql 


( 1 . 11 ) 


The steady-state frequency response at the junction x = 0 of fig. 1 now becomes 


T(0) = 


1 

G(co) 


(1 -sech ql). 


( 1 . 12 ) 


The term 1 - sech (ql) in Eq. (1.12) corresponds to the steady-state 
temperature distribution for a constant ambient fluid temperature Tf > T c 
derived by Scadron and Warshawsky (1952). It is also interesting to note that the 
first order response can be derived from Eq. (1.12) in the limit of large wire length 
1 or 


limT(O) = — - — 
1_> “ G(co) 


(1.13) 


where Eq. (1.13) represents the frequency response of a uniform thermocouple 
wire with no axial heat conduction. 
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The steady-state frequency at the thermocouple junction x = 0 is normally 
characterized graphically in the form 

T(0) = IT(0)le* (l- 14 ) 

where I T(0) I is the amplitude ratio and <j> is the phase angle, in the latter case, 


the phase angle in degrees is 


<J> = 57.3 tan 


-1 


ImT(O) 
Re T (0) 


(1.15) 


where Im[T(0)] and Re[T(0)] are the imaginary and real parts of T(0), respectively. 


1.2.2 Previous Work 

Scadron and Warshawsky (1952) in their pioneering work considered the 
frequency response of a uniform thermocouple wire as shown in fig. 1. 
Introducations for the dependent variable T w in Eq. (1.1) as suggested by Carslaw 
(1945), Scadron and Warshawsky obtained expressions for the steady-state 


frequency response in the form of the amplitude ratio 


IT(0)l = 



( Y 

s? + 

-1 



1+ ^ 
l 0)„ 


(1.16) 


and phase angle 


—i/ toSo . , ® \ 

<J> tan ( — -*-) - tan ( — ). 
WpSt co n 


(1.17) 


Here, the infinite series Sj and S 2 that appear in Eqs. (1.16) and (1.17) are 
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s, = i-£. 


i=i 


r 

1 

dj(dj-l) 

d ' + 

f— T 

) 

l w J . 


and 


s 2 =£ 


c, 


j-1 


d: “ 1 


d^ + 


^ CO ^ 

V®ny 


where the constant that appear in the terms of Si and S 2 are 
c i = ^-[1 - cos(jn)]sin p K ^ X ^ j 


and 


d,=i + Mi. 

i 41 2 


(1.18a) 


(1.18b) 


(1.19a) 


(1.19b) 


Properties of these solutions derived by Scadron and Warshawsky are discussed 
in Sec. III. 


1.2.3 Non-Uniform Wire 

If the wire diameters or the physical properties such as thermoconductivity 
are different across the thermocouple junction as shown in fig.2, Eq.(1.6) must be 
solved in the separate regions 1 and 2 on both sides of the junction. This 
solution is subject to the boundary conditions of continuous temperature and 
equal rates of conductive heat transfer at the junction x = 0 or 

= T 2 (l ) = 0 (1.20a) 

71(0) = 7 2 (0) (1.20b) 
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(1.20c) 


KD ]«m = ^ 

dx 


2 ^(0) 

dx 


In region 1, the general solution to Eq. (1.6) for the frequency response can be 
written in the form 

T x (x) = A x sinh q x (l + x) + B x cosh + + ’ (1 ' 21) 

Since T x (-/) = 0 from boundary condition Eq. (1 .20a), one finds for the constant 
Bt 


0 = B, + 


G^co) 


Substituting for Bi into Eq. (1.21), one obtains 

TAx) = A sinh< 7 ,(/ + x) + — — -[l — cosh (ft (/ + *)] (1.22) 

1 1 G x {co) 

where the parameters in Eq. (1.22) reflect the material properties and wire 
diameter of region 1 or 


= 


P\ C )D\ 


4 ^-, G x (co) = \ + i 


co, 


CD 

— [ Q\ 

V 


1 


a, (a) 


r, 


Similarly in region 2, the frequency response can be shown to be of the form 

T 2 (x) = A 2 sinh q 2 (l-x) + — — -[l-cosh q 2 (l-*)] (1.23) 

C/2 yu*) 

where the constants Ai and A 2 in Eqs. (1.22) and (1.23) above must be determined 
by the simultaneous solution to the boundary conditions Eqs. (1.20b) and (1.20c). 
Substituting Eqs. (1.22) and (1.23) into boundary condition (1.20b), one 

obtains a linear equation for the constants of integration Ai, A 2 or 

A x sinh q x l - A 2 sinh q 2 l = ■ (1-24) 
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Similarly, substituting Eqs. (1,22) and (1.23) into boundary condition (1.20c), 
obtains a second linear equation for A j , A 2 or 


one 


A Q c °sh qj + A 2 cosh q 2 l = R 2 . 
Here, the complex constant 

Q = M,V,)/(* 2 D 2 V 2 ) 

and the complex constants R,, R 2 are equal to 

r> __ ( ^ 1 v 1 | 

' ” A ~ (V ^ COsh 9,1-— cosh qj 

and 


(1.25) 


(1. 26) 


(1.27) 


= ^si„h„,, + J-si„h ?2 l. (,.28) 


Solving Eqs. (1.24) and (1.25) for A, and A 2 , one obtains the determinate 


system 


A = 


A - sinh< 7 2 / 
R 2 cosh q 2 l\ 


\ = 


DET 

sinh( 7 ,/ /?, 

QcoshqJ R ' 


DET 


(1.29) 


(1.30) 


where the determinate in the denominator is equal to 

DET = sinh qi / cosh q 2 / + Q cosh q] / and q 2 / 

Thus, the constants A : and A 2 in Eqs. (1.22) and (1.23) become 

A = -h f + tanh *,) 

smh qj + Q cosh qj tanh q 2 l 

and 


(1.31) 


.(1.32) 
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(1.33) 


Az — 


Hi^ftanh q,l - tanh & ) + tanh ^ + Q' ™ h 

Gi \ . — . ^ — — 


sinh q 2 l + Q' cosh q 2 l tanh qj 


where the complex constant in A 2 is 

Q' = 1 /( 2 - 

Finally, the steady-state frequency response at the thermocouple junction x = 0 
can be determined by substituting the value for Ai above or Eq. (1.32) into the 
general solution for T, (0) given by Eq. (1.22). One obtains after simplification 

Q^3l { , - coth «,/ tanh &)+■ ' coth qj tanh ^ ^ 

C/| V ^ 

T i(°) = Q coth q x l + coth q 2 l 

where it can be shown that T,(0) = T 2 (0) at the thermocouple junction. 

It is interesting to note that one can determine a frequency response for a 
non-uniform wire of infinite length from Eq. (1.34) in a manner similar to the 
derivation of Eq. (1.13). Taking the limit of infinite wire length in T,(0> above, 
one obtains the simple expression 


Q 1 

— + — 

G, G 2 

limT,(°) = 


(1.35) 


Equation (1.35) represents the frequency response for a non-umform wire which 
isl. analogous to the well known first order response given by Eq. (1.13) for a 

uniform wire of infinite length. 
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1.3 RESULTS 

The amplitude ratio and phase angle of the thermocouple frequency 
response were plotted graphically for the case of a uniform wire with the 
geometry shown in fig. 1. In this case, average properties of a type B or 
Pi / 6% Rh-Pt / 30% Rh were used since the material properties were nearly equal 
across the thermocouple junction. The wire dimensions, properties and gas 
conditions are listed in table KTouloukian et al., 1970). 

The amplitude ratio and phase angle were also plotted for a composite or 
non-uniform thermocouple wire with the geometry shown in fig. 2. In this case, 
different wire diameters and material properties were used including a type B 
thermocouple described in Table 1 in addition to a type T or copper - constantan 
described in table 2 (Touloukian et al., 1970). 

The form of the convective heat transfer coefficient h that appears in the 
computation of the natural frequency (O n defined in Eq. (1.8) was determined 
from the expression (Scadron and Warshawsky 1952) 

Nu= .431 Re^ (1.36) 

where Nu (= hD / k f ) \ s the Nusselt number, k f is the thermal conductivity of the 
ambient fluid and Re(=vD/u / ) is the Reynolds number of the thermocouple 
wire. Here, v and v f are the fluid velocity and kinematic viscosity, respectively. 

It should be noted that the convective heat transfer coefficient h a D xn and the 
natural frequency of a thermocouple wire for given material properties 
co n a D ~ V1 . 
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1.3.1 Uniform Wire 


The amplitude ratio |T(0)| at the wire junction for the steady-state frequency 
response Eq. (1.12) is computed in fig. 3. This assumes a type B thermocouple 
wire with the dimensions listed in table 1. In this case, the average material 
properties listed in table 1 were used since the theory assumes that the 

thermocouple is uniform across the junction. In fig. 3 the amplitude ratio |T(0)| 

-4 .54 for low frequencies o) _> 0. This lower limit corresponds to the attenuation 
expected for a thermocouple wire exposed to a steady ambient gas temperature Tf 
> T 0 as derived by Scadron and Warshawsky (1952). 

The amplitude ratio Eq. (1.16) predicted by Scadron and Warshawsky (1952) 
is also graphed in fig. 3. It was found that the series Eqs. (1.18a) and (1.18b) 
converged slowly such that the final result was in error by ~ 11% at low 
frequencies after retaining the first 10 terms. The first 50 terms of the series 
expressions for Siand S 2 reduced the error to 2.3% while retaining 100 terms as 
was done in fig. 3 limited the error to 1%. The final result is identical, however, 
to the amplitude ratio predicted by Eq. (1.12) and derived in the present chapter. 

The amplitude ratio for a type B thermocouple wire of infinite length from 
the first order response Eq. (1.13) is also indicated in fig. 3. As expected, the first 
order response or Eq. (1.13) asymptotically approaches the amplitude ratio 
predicted for a wire of finite length in the limit of large frequencies (O I co n » 1 
where the time scale of the ambient temperature fluctuation is too small for axial 
heat conduction to occur. 
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The phase angle <j> for all three theories is computed in fig. 4. The phase 
angle varies over the range 0 > O > - nj2 and approaches the lower limit of 
- nil for large frequencies co / co n » 1 . it is apparent that the effect of axial heat 
conduction is to reduce the magnitude of the phase angle. 

The spacial variation of the amplitude ratio |T(*)| derived from Eq. (1.11) is 

graphed in fig. 5. The present theory approaches the correct limit of |T(±/)| = 0. 
However, the theory of Scadron and Warshawsky oscillates near the end of the 
w ’ re x — » ±/. In the latter case the amplitude ratio equals the value 

l/Jl+(co/co n ) for x = ± l since Si = 1 and S 2 = 0 in Eqs. (1.18a) and (1.18b), 
respectively. An incorrect value for the amplitude ratio near the ends of the 
wire would not ordinarily present a problem unless one wished to adapt the 
theory to the common configuration of a supported thermocouple wire with 
heat transfer down the support legs (Forney and Fralick, 1991). 

1.3.2 Non-Uniform Wire 

The amplitude ratio |7^ (0)| from the steady-state frequency response Eq. 
(1.34) is graphed in fig. 6. This plot assumes a type T thermocouple wire with the 
dimensions listed in table 2 . An important parameter that affects the relative 
magnitude of the terms in Eq. (1.34) and thus the relative importance of the 
thermocouple wire segments in regions 1 or 2 in the schematic of fig. 2 is the 
complex dimensionless ratio Q = ktfqj k 2 D 2 2 q 2 . We note here that Q assumes 
different values in the limit of either low or high frequencies. Given that the 
Nusselt number is related to the wire Reynolds number by Eq. (1.36), one obtains 
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for low frequencies 


Q{co = 0) = Qo = 

while for high frequencies 

Q(o) = °°) = (L = It- 


.5 / 

‘ ( ' 
\ D 2 J 


1.25 


(1.37) 


(V 

.5 / \2 

(£l) 

UJ 

U> 2 J 


(1.38) 


The amplitude ratio |r t (0)|-> .228 for low frequencies w-^O in fig- 6. The 
lower limit corresponds to the attenuation predicted by Scadron and 
Warshawsky (1952) for a non-uniform thermocouple wire as shown in fig. 2 that 
is exposed to a steady-state ambient gas temperature T f > T 0 . The real expression 
for the temperature distribution developed by Scadron and Warshawsky for this 


case is 


TMzIa = 


(1.39) 


where 


¥ = 


, ly* , ly 

sech cosh — - 

2 2 - 


ly~ . ly* 
1- m^tanh-f- coth— - 
JL £ 


1 - n\ 2 sinh ly csch ly* 
Here, the parameters y + and y" are defined as 


r+ = T = + T = ’ 


y~ = 




and 


m n = 


0,-1 

Qo + 1 


(1.40) 


(1.41a) 


(1.41b) 
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The amplitude ratio fc(0)| from Eq. (1.35) for a type T thermocouple wire 
with the same wire dimensions D 2 / D, = 2 as the solid curve in fig. 6 but of 
infinite length / -» „ is also computed. As indicated, the amplitude ratio at low 
frequencies ru -> 0 for a non-uniform thermocouple wire of infinite length 
approaches fc(0)|->l. This result is identical to the first order frequency response 
Eq. (1.13) for a uniform wire of infinite length as is shown in fig. 3. For the non- 
uniform wire, however, axial heat transfer does occur at the junction because of 
dissimilar material properties. At large frequencies co/o>,» 1, however, the 
amplitude ratio derived from Eq. (1.35) for the wire of infinite length 
asymptotically approaches the value predicted by Eq. (1.34) for a wire of finite 
length as shown in fig. 6. The latter result is similar to the features of fig. 3 at 
large frequencies for the case of a uniform thermocouple wire. These latter 
features reflect the fact that the time scales of the ambient gas temperature 
fluctuations at high frequency are too small for axial heat conduction to occur. 

The phase angle <j, for both expressions Eqs. (1.34) and (1.35) is graphed in 

fig- 7. The phase angle varies over the range 0 >^>-n /2 and approaches the 
lower limit of - n/2 for large frequencies . It is apparent that the effect of axial 
heat conduction is to reduce the magnitude of the phase angle. 

The amplitude ratio |7j(0)| derived from Eq. (1.34) for a type B thermocouple 
wire is computed in fig. 8. Since the material properties in this case are nearly 
equal on both sides of the junction, the indicated difference in the amplitude 
ratio is the result of different wire diameters across the junction. These results 
are strongly influenced by the magnitude of the factor Q„ defined by Eq. (1.37). If 
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the parameter Qo - 1 as is the case for the upper curve in fig. 8, both regions 1 and 
2 in the schematic of fig. 2 are equally weighted in the expression for the 
frequency response given by Eq. (1.34). However, increasing the wire diameter D 2 
in region 2 means Qo is reduced in magnitude such that Qo < 1- Now the 
amplitude ratio is more strongly influenced by the properties of the larger wire 
in region 2. The phase angle is plotted in fig. 9 for the same wire configuration. 

The amplitude ratio ft Ml is also computed as a function of the axial 
distance down the type B thermocouple wire in fig 10. This was accomplished 
by substituting the expression for the constant A, given by Eq. (1.32) into the 
general solution for the frequency response Eq. (1.22). As discussed above, if the 
ratio Qo < 1 because of an increase in D 2 , a large increase in the rate of heat 
transfer occurs in region 2 down the larger wire. This additional axial heat 
conduction reduces the amplitude ratio at the thermocouple junction compared 
to the case of Qo - 1 for the wire geometry represented by the upper curve in fig. 

10 . 

The amplitude ratio derived from Eq. (1.34) for a type T thermocouple wire 
is graphed in fig. 11. Since the thermoconductivity of copper in region 1 is more 
than an order of magnitude larger than the thermoconductivity of constant in 
region 2, the parameter Qo > 1 for D, /Dj - 1. Increasing the diameter of region 2, 
however, has a small effect on the amplitude ratio since Qo > 1, which remains 
unchanged, and the properties of region 1 continue to dominate the rate of heat 
transfer. The phase angle for the same configuration is computed in fig 12. 
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The amplitude ratio ft(x)| is also computed as a function of the axial 
distance down the type T thermocouple wire in fig. 13. As discussed above, Qo> 
1 for both wire configurations. Thus, the large rate of axial heat transfer down 
the copper wire in region 1 gives nearly equal spacial temperature profiles for 
both wire configurations and this fixes the junction temperature. 
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1.4 NOMENCLATURE 


A = 
B 

c = 

c i = 
dj = 
D = 
G = 
h 

i = 
k = 
k f = 
/ 

mu = 
Nu = 

Q = 

Qo = 

Q- = 

q = 

Re = 
Si = 


constant of integration 
constant of integration 

material specific heat (J~S m 1 - K 
term in infinite series Si , S 2 
term in infinite series Si , S 2 
thermocouple wire diameter (cm) 

1 +/(©/©.) 

heat transfer coefficient [J -cm 2 - s 


unit imaginary number 

material thermoconductivity {J -cm 1 

gas thermoconductivity ( J — cm -s 
length of thermocouple wire (cm) 

(&-!)/(& + !) 



Nusselt number (= hD/kf) 
(k x D x \ )!(k 2 D 2 q 2 ) 


{k i /k 2 ) 5 {D l /D 2 f 75 
(*, / k 2 )\D, / D 2 f 

(G/r ) 5 

Reynolds number (= vD/ v f ) 
infinite sum 
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infinite sum 


time(s) 

frequency response 

amplitude of periodic gas temperature (°^) 

gas temperature (* ) 

steady-state gas temperature (°^) 

complex amplitude of periodic wire temperature(°^) 

local wire temperature (°A') 

gas velocity (cm - s ~' ) 

axial distance from center of wire (cm) 

Greek Symbols 

thermal diffusivity (cm 2 -s~ x } 

oc / C0 n (cm 2 } 



kinematic viscosity of gas (cm 2 

angular frequency (s -1 ) 

natural frequency of wire (=4 h / pcD)(s~'} 
phase angle (degrees) 


normalized temperature at junction 

-3 

material density (gm - cm ) 
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1.6 TABLES AND FIGURES 
Table 1 - Uniform Wire Conditions 


Dimens ions (cml 
D = .0076 


Air Properties 
T 0 = 300 °K 


l = 


Pt - 6% Rh 
Pt - 30%Rh 
Average 


p = 1 atm 

v = 50 m/sec 
Propertie s of Type B 


I \ 

f cm 2 ^ 

J — 

a\ 

{cm'-°K) 

Vsec 

2.73 

.238 

2.86 

.190 

2.8 

.214 


Pj (cml 

.0076 

.0076 


Table 2 - Properties of Non-Uniform Wire 


Dimensions (cm) 
D ?(cm) D 2 /Pj /( cr o) 
.0076 1 1 

.0152 2 -1 


Wirp Location 

Region Type B T ype T 

1 Pt-6%Rh Copper 

2 Pt-30%Rh Constantan 


Properti es of Type T 




a 

cm 
sec _ 

Copper 

3.44 


1.16 

Constantan 

3.48 


.067 
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Fluid Flow 



I - / i 77 

Fig. 1 Schematic of uniform wire 


Fluid Flow 



Fig. 2 Schematic of non-uniform wire 
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Fig. 5 Amplitude ratio vs axial distance for uniform wire. 
Dotted line, Eq. (11). 
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Fig. 10 Amplitude ratio vs axial distance foi 
material properties. Each curve is de; 
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2.0 SUPPORTED THERMOCOUPLE WIRE: THEORY AND EXPERIMENT 

Theoretical expressions are derived for the steady-state frequency response 
of a supported thermocouple wire. In particular, the effects of axial heat 
conduction are demonstrated for both a supported one material wire and a two 
material wire with unequal material properties across the junction. For the case of 
a one material supported wire, an exact solution is derived which compares 
favorably with an approximate expression that only matches temperatures at the 
support junction. Moreover, for the case of a two material supported wire, an 
analytical expression is derived that closely correlates numerical results. 

Experimental data were taken with a type K supported thermocouple. The 
test thermocouple was constructed with dimensions to demonstrate the effects of 
axial heat conduction assuming constant physical properties across the junction. 

2.1 INTRODUCTION 

In the present chapter, the theoretical steady-state frequency response of a 
supported thermocouple wire has been calculated to include the effects of axial 
heat conduction. These solutions, which represent an extension of earlier work 
(Forney and Fralick, 1991) are derived for both a supported thermocouple wire 
with equal physical properties across the junction (e.g., roughly the same 
thermoconductivity, etc.) and a supported wire with unequal properties across the 
junction. Solutions are presented in the form of the amplitude ratio and phase 

angle for both cases. 
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2.2 THEORY FOR ONE MATERIAL THERMOCOUPLE 

The steady-state frequency response of a thermocouple wire will be 
developed with the following assumptions: (a) the amplitude of the fluctuating 
fluid temperature is small relative to the mean absolute temperature (b) the 
thermocouple dimensions are small relative to the size of the turbulent eddies or 
enclosure dimensions (c) radial temperature gradients in a wire cross section can 

be neglected and (d) radiative heat transfer can be neglected relative to conduction 
and convection. 

In this section the geometry of Fig. 1 is considered where the material 
properties of thermal conductivity k, specific heat c and wire density pare assumed 
to be equal on both sides of the thermocouple junction. If the probe is immersed 
in a flowing fluid, the expression for the local conservation of energy in the 
thermocouple wire becomes (Scadron and Warshawsky, 1952) 

dt dx 2 + P cd( T *~ T w ) ( 2 .d 

where a = k/pc is the thermal diffusivity of the wire, T g is the ambient fluid 
temperature, h is the convective heat transfer coefficient, D is the wire diameter 

and T w is the local wire temperature measured along the axis at a distance x from 
the centerline (Fig. 1.) 

The wire and fluid temperatures are measured relative to the mean fluid 
temperature T 0 . The ambient fluid temperature is taken to be a mean temperature 
together with a sinusoidal varying deviation from the mean, 

T g (t) = T 0 + T f e'“< (2.2) 

where ® is the angular frequency of the ambient temperature. Since Eq. (2.1) is 


39 



linear, we now seek a solution for the local wire temperature of the form 
(Hildebrand, 1976) 

T. = T 0 + 7 i(x)e-. C ' 3) 

Referencing all temperatures with respect to the mean gas temperature T 0 and 
normalizing with respect to the amplitude of the fluctuating ambient fluid 
temperature T f , one defines a local normalized steady-state frequency response 

f(x) for the thermocouple wire of the form 


ZWo = = T(xy ml . 


(2.4) 


T, 


*/ i 

Substituting Eqs. (2,2) and (2.3) into Eq. (2.1), one obtains an ordinary 
differential equation of the form 


icoT = a + 
dx 


(2.5) 


where T = T(x) is the frequency response and the vector notation will be dropped 
for simplicity. Thus, for the geometry of fig. 1, one seeks a solution to the non- 
homogeneous linear second order differential equation for the dependent variable 

T of the form (Forney and Fralick, 1990) 

y T" - G(co)T = - 1. (2 - 6) 

The general solution to Eq.(2.6) can be written in the form (Hildebrand, 1976) 

1 


T(x) = A sinh qx + Bcosh + 


C(o) 


(2.7) 


where the parameters in Eqs. (2.5) and (2.6) are defined as 


m - — r = — , G(<») = 1 + i — 

“■ pcD' 7 (O. I®., 


( 2 . 8 ) 
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while in Eq. (2.7) the constants A and B are complex, 1/G(<y) represents the 
particular solution and the parameter 



2.2.1 Approximate Solution 


Assuming that the material properties are constant across the junction and 
that the wire diameters are D] and D 2 in regions 1 and 2, respectively, Eq. (2.6) is 
subject to the boundary conditions 


T x (l) = T 2 (l) = T a (2.9 a ) 

T 2 (l + L) = 0 . (2.9b) 

In this case, we seek a simple approximate solution that neglects the heat transfer 
at the interface between regions 1 and 2 at x = ± l where the parameters in Eqs. 
(2.5) and (2.6) are defined in terms of the wire diameters in each region. A similar 
approach will be used in a later section to obtain an approximate solution for the 
case in which the material properties of the two elements of the thermocouple are 
distinctly different. Hence, in region 1 


a, = 


_ 4/z, 

pcD, 


y - — > g, (oj) = i+ /(— ) 


where ty, is the natural frequency of the wire in region 1 of fig. 1. 

The solution to Eq. (2.6) for the one material wire on both sides of the 


junction in region 1 can be written in the form 


7JU) = A ] sinh q t x + B l cosh qp + 


1 

G,m ‘ 


( 2 . 11 ) 
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Substituting x = ±l in Eq. (2.11), the boundary conditions in Eq. (2.9) yield values 
for the constants 


cosh qj 


( 2 . 12 ) 


Thus, one obtains a steady-state temperature distribution for the wire in region 1 


of Fig. 1 in the form 


TM = 7T 1- 


1 f cosh<7,x 


G ] l cosh^/ 


cosh^jc 

cosh( 7 ,/ 


(2.13) 


We now seek a solution in region 2 that satisfies the boundary conditions of 
Eq. (2.9). Since the temperature is symmetric about x = 0, it is convenient to define 
a continuous steady-state temperature distribution for the large wire of diameter 

E >2 over the entire region -(l+L) < x < {l+L) or 


TJx) = 


cosh < 7 2 x 
cosh^ 2 (/ + L) 


(2.14) 


Since T a = T 2 (l), one obtains the boundary value from Eq. (2.14) 


1 L cosh </ 2 / 
a G 2 l coshq 2 (l + L) 


(2.15) 


Substituting the value for T a aix-l of Eq. (2.15) into Eq. (2.13), the 
approximate temperature distribution in region 1 becomes 


1 ( cosh<7,x^| ( l _ cosh <y 2 / ^j cosh^x 

Ti(X) = a l 1 ” cosh qj ) g 2 ^ cosh^ 2 (/ + L) J cosh <7,/ 


(2.16) 


Thus, the approximate frequency response at the thermocouple junction (x = 0) for 


the one material wire becomes 
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(2.17) 


1 

cosh <7/ 

The steady-state frequency at the thermocouple junction x=0 is normally 
characterized graphically in the form 

7(0) = |r(O)|e i0 (2.18) 

where |7(0)| is the amplitude ratio and <j> is the phase angle. In the latter case, the 
phase angle in degrees is 


^ <0) = i 


f 


1 - 


1 _ 

cosh^,/ 




1 - 


i 1 ) 


cosh<7 2 / 

cosh (^ + L) 


O = 57.3 tan' 1 


Im 7(0) 
Re 7(0) 


(2.19) 


where Im[T(0)] and Re[T(0)] are the imaginary and real parts of T(0), respectively. 


2.2.2 Exact Solution 

If the boundary conditions listed in Eq.(2.9) include equal rates of conductive 
heat transfer at the interface between the thermocouple and support wires at 
x = ± /, the exact solution is subject to 


II 

f-r 

(2.20a) 

1 dx 2 dx 

(2.20b) 

7 2 (/+L) = 0 . 

(2.20c) 

Since the solution to Eq. (2.6) in region 1 is of the form 


T x (x) = A l sinh< 7 ,jc + fijcosh^x + — 

G, 

(2.21) 


where by symmetry 7^ (/) = 7^ (-/), one obtains A i = 0. Thus, the form of the 
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solution in region 1 becomes 


7^(x) = Z^coshc^x + — . 


( 2 . 22 ) 


For region 2, where the spatial coordinate is in the range l<x<(l + L ), it is 
convenient to write 


T 2 (x) = A 2 sinh q 2 (/ + L - x) + B 2 cosh<? 2 (/ + L — x) + — 


(2.23) 


From the boundary condition T 2 (l + L) = 0, one obtains B 2 = -1 /G 2 or the form of 
the solution in region 2 becomes 


T 2 {x) = A 2 sinh< 7 2 (/ + L — x) + — [l - cosh< 7 2 (/ + L - *)] . 

G 2 


(2.24) 


Substituting Eqs. (2.22) and (2.23) into boundary condition (2.20a), one obtains 
a linear equation for the constants B,, A 2 or 

B, cosh <7,/ - A 2 smhq 2 L = -J-[l - cosh<? 2 L] - -J- . 

C/2 C/j 

Similarly, substituting Eqs. (2.22) and (2.23) into boundary condition (2.20b), one 
obtains 


a second linear equation for B, , A 2 or 


1 . 


B x Qs\nhq x l + A 2 coshq 2 L = — sinh< 7 2 L 


(2.26) 


where the complex constant Q is defined as 


Q^. 

D\q 2 


(2.27) 


Solving Eqs. (2.25) and (2.26) for B] and A 2 , one obtains the determinate 
system 
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5 , = 


-fl -cosh<7 2 L]- 


1 


C, 


sinh<7 2 L 


DET 


— sinh< 7 2 L| 
cosh q 2 L 


(2.28) 


and 


^2 — 


cosh < 7,1 
QsinhqJ 


^-[l-cosh<7 2 L]-^ 


1 


G i 


sinh< 7 2 L 


DET (2.29) 

where the determinate in the denominator is equal to 

DET = coshq,/ cosh</ 2 L + £ 2 sinh< 7 ,/sinh< 7 2 L. (2.30) 

Thus, solving for the constants B] and A 2 from Eqs. (2.28) and (2.29) and 
substituting B] into Eq. (2.22), one obtains an exact expression for the steady-state 
temperature distribution in the form 


T x (x) = 


— (cosh< 7 2 L - 1 ) - — cosh< 7 2 L 
G 2 G, 


cosh < 7 ,* 


1 

+ — . 

G, 


cosh < 7 ,/ cosh< 7 2 L + Qsinh^,/ sinh< 7 2 L G, (2.31) 

Thus, the steady-state frequency response at the thermocouple junction x=0 
becomes 


— (cosh< 7 2 L- 1 ) — cosh< 7 2 L 


T( 0 ) = 


G, 


1 


cosh < 7 ,/ cosh< 7 2 L + < 2 sinh< 7 ,/ sinh< 7 2 L G, 


(2.32) 
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2.3 THEORY FOR TWO MATERIAL THERMOCOUPLE 

Certain types of thermocouples have distinctly different material properties 
across the junction. For example, a copper-constantan thermocouple has a 
thermal conductivity on the copper side that is more than an order of magnitude 
larger than constantan. In this case, the expressions developed in the previous 
section for the frequency response are in error since unequal material properties 
would provide an asymmetric temperature profile. 

Referring to fig. 2, the thermocouple schematic now has four distinct regions 
that are distinguished by either different wire diameters or physical properties. For 
example, on the left side of the schematic of fig. 1 the thermal conductivity, density 
and specific heat have the values ki, p, and ci , respectively, while on the right side 
of the schematic the material properties are k 2 , p 2 and C 2 - 


2.3.1 Temperature Distribution for Small Wire 

Since the differential equation describing the steady-state frequency response 
Eq. (2.6) applies to all regions of the schematic of fig. 1, the steady-state frequency 
response in regions 1 and 2 are given, respectively, by the expressions 


T^x) = i4,sinh<7jX + B ] cosh q { x + — 

Gj 

T 2 (x) =A;,sinh^ 2 x + B 2 cosh q 2 x + — . 

G 2 

Here, the four constants designated by A and B in Eqs. 
determined by the four boundary conditions 


(2.33) 

(2.34) 

(2.33) and (2.34) are 
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7i(0) = T 2 ( 0) 


(2.35a) 


k]D >*m = kiD M 


dx dx 

n-D = t, 
m= T h . 


(2.35b) 

(2.35c) 

(2.35d) 


Solving for the four values of the constants designated by A and B in Eqs. 
(2.33) and (2.34), one obtains an expression for the steady-state temperature 
distribution in region 1 of the form 


7'iU)"4- + l sinh fl F i(- r+/ ) 

C/, 


( \ 1^ A 

V G 2 G, y 


cosh q 2 l + 


( 


+ 


A 


T - — 

v* G „ 


(2 r cosh sinh - sinh qx cosh ql 



(2.36a) 


and in region 2 


T 2 (x) = + sinh <7 2 (/_x) l 

2 


vS c „ 


cosh^/ - 


V - G 

« G 

v ly 


V-G 

b G 
V 2J 


[cosh q_x sinh ql + ^ cosh sinh — 


(2.36b) 


Here, it should be noted that the two functions above can be obtained from each 
other by interchanging the subscripts 1 and 2, the values a and b and by changing 
the sign of x. Also, in Eqs. (2.36a) and (2.36b) the parameters 

A; = <2; cosh q x l sinh q 2 l + cosh q 2 l sinh qj 


and 


a = 


* I <?1 

k 2 q 2 
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2.3.2 Temperature Distribution for Large Wire 

In this section a solution is sought for the large wire in regions 3 and 4 of the 
schematic of fig. 2. To simplify the analysis, a continuous solution is derived for 
the steady-state frequency response over the entire region -(/ + L) < x < (/ + L). This 
solution must satisfy the boundary conditions 


7(0) = 7(0) 

(2.37a) 

dT(0) dT( 0) 

kD 2 3 = k D 3 4 

(2.37b) 

1 3 dx 2 * dx 

I 
1 

II 

o 

(2.37c) 

+ 

II 

o 

(2.37d) 


In principle, it is now possible to derive a continuous solution for T 3 (x) over 
the range — (/ + L) < x < 0 and for T 4 (x) over the range 0 <*<(/ + L). These solutions 
are obtained from Eqs. (2.36a) and (2.36b) by replacing the subscripts 1 — » 3 and 2 — » 4 
and by redefining the parameters / / + l and setting the boundary conditions 

T a = T b = 0. Moreover, the boundary condition T a =T 3 (-l) is obtained from the 
resulting expression for T 3 (x) by substituting x = -l while the boundary condition 
Tb = T 4 (l). Thus, one obtains values for T a , Tb that appear in Eqs. (2.36a) and (2.36b) 
in the form 

T a = -5- + {sinh f-J — -J-lcosh^ 4 (/ + L)--J- 

L 3 Cr 4 C / 3 J G 4 

- —\Q, cosh q 3 l sinh q 4 (l + L) + cosh# 4 (/ + L)sinh<7 3 /l}f— 

C 3 *\A J ( 2 . 38 ) 


48 


and 



{Q s sinh^L 


Y i 

l 3 

_V^3 

cj 


cosh< 7 3 (/ + L)-- 


[cosh qj sinh q 3 (l + L) + Q s cosh q 3 {l + L)sinh qj]\ 


f i A 


V^2 ) 


(2.39) 


Here, the parameters are defined as 

A 2 = Q s cosh q 3 (/ + L) sinh q 4 (l + L) + cosh q 4 (l + L) sinh q 3 (l + L) 


and 


Qs = 


Mi 

k 2 q 4 


2.3.3 Frequency Response 

The steady-state frequency response for the two material supported 
thermocouple is obtained from Eq. (2.36a) or (2.36b) by setting x = 0: 


r(0) = — \q, sinh q 2 l 


T a + yr{c osh<7,/-l) 

Cr, 


+ sinh^J 


T b +-^r{coshq 2 l-l) 

CTt 


(2.40) 


where T a and Tb are given by Eqs. (2.38) and (2.39). 

It should be noted that the steady-state frequency response Eq. (2.40) provides 
a wire temperature that is continuous everywhere and conserves the heat flux at 
the junction x = 0. This represents an approximate solution since the heat flux at 
the interface between the large and small wires x = ±l has been neglected. 


2.4 RESULTS 

The amplitude ratio and phase angle of the thermocouple frequency response 
were plotted graphically for the case of a one material wire as shown in the 
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schematic of fig. 1. In this case, average properties of a type B or 
Pt / 6% Rh- Pt / 30% Rh were used since the material properties were nearly equal 
across the thermocouple junction. The wire dimensions, properties and gas 
conditions are listed in table l(Touloukian et al., 1970). 

The amplitude ratio and phase angle were also plotted for a two material 
thermocouple wire as shown in the schematic of fig. 2. In this case, a type B 
thermocouple described in table 1 was used in addition to a type T or copper - 
constantan described in table 2 (Touloukian et al., 1970). 

The form of the convective heat transfer coefficient h that appears in the 
computation of the natural frequency co n defined in Eq. (2.8) was determined from 
the expression (Scadron and Warshawsky 1952) 

Nu = .485 Re^ Pr^ (2.41) 

where Nu(= hD / k f ) is the Nusselt number, kf is the thermal conductivity of the 

ambient fluid, Pr(= u 7 / a) j s th e Prandl number, and Re(=vD/i> / ) is the 
Reynolds number of the thermocouple wire. Here, v and v / are the fluid velocity 
and kinematic viscosity, respectively. It should be noted that the convective heat 
transfer coefficient h a D 1/2 and the natural frequency of a thermocouple wire for 
given material properties (o n a D _3/2 . 

2.4.1 One Material Thermocouple 

The amplitude ratio |T(0)| at the wire junction for the steady-state frequency 
response derived from Eqs. (2.17) and (2.32) is shown in fig. 3. This assumes a type 
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B thermocouple wire with the dimensions listed in table 1. In this case, the 
average material properties listed in table 1 were used since the one material 
theory assumes that the properties of the thermocouple wire are uniform across 
the junction. It is evident in fig. 3 that the amplitude ratio derived from the 
approximate expression Eq. (2.17) is nearly identical to the exact derivation Eq. 

(2.32) . Thus, it appears that the conservation of heat flux at the interface between 
the small and large wires of the schematic of fig. 1 is of secondary importance. 

Also included in fig. (3) are the numerical computations of Stocks (1986). 
These solutions represent explicit finite difference solutions to the one 
dimensional unsteady heat transfer equation as shown in Eq. (2.1). The small 
deviation of the numerical results from the exact solution at low frequency in fig. 
(3) is apparently due to the unsteady character of the numerical results. Similar 
computations of the phase angle <j> for the type B thermocouple are represented in 
fig. 4. As indicated, the phase angle varies over the range 0>0>-7t/2 and 
approaches the lower limit of -n/2 for large frequencies » 1 • 

The spacial variation of the amplitude ratio |T(x)| derived from Eqs. (2.17) and 

(2.32) is graphed in fig. 5. These computations were made at an angular frequency 
of (o/co n = 0.1 for the type B thermocouple. As evident in fig. 5, the difference 
between the exact and approximate expression is somewhat exaggerated at a very 
low frequency. Nevertheless, the error represented by the approximate solution is 
less than 7% over the length of the thermocouple. As stated earlier, matching the 
heat flux at the interface between the small and large wire at x = ± / in the 
schematic of fig. 1 appears to be of secondary importance in relation to providing a 
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continuous temperature profile along the wire. 

The amplitude ratio |T(0)|, phase angle cj> and spacial variation p (*)| are also 
plotted in figs. (6), (7) and (8), respectively, from the steady-state frequency response 
represented by Eqs. (2.17) and (2.32). These results represent a type B thermocouple 
with a smaller diameter ratio D 2 jD x = 2 (see table 1 for dimensions). As indicated 
in figs. (6), (7) and (8), the approximate and exact solutions represent comparable 
results in all three cases. Thus, the diameter ratio of the large and small wire 
appears to have little effect in the favorable comparison between the approximate 
and exact solutions representing the one material steady-state frequency response. 

2.4.2 Two Material Thermocouple 

The amplitude ratio (0)| at the wire junction for the steady-state frequency 
response derived from Eqs. (2.17) and (2.40) is shown in fig. 9. This assumes a type 
B thermocouple wire with the dimensions listed in table 2. In this case, the 
average material properties listed in Table 2 were used for the amplitude of the 
frequency response derived from the one material solution of Eq. (2.17). Also 
plotted in fig. 9 is the amplitude ratio derived from the two material solution of 
Eq. (2.40). In the latter case, the individual material properties listed in Table 2 
were used. 

As expected, the amplitude ratio for the steady-state frequency response of a 
type B thermocouple is nearly identical with either the one material or two 
material approximate solutions. This is a consequence of roughly equal material 
properties across the junction for type B thermocouples. This plot also validates 
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the two material approximate solution Eq. (2.40) and the values of the boundary 
conditions for T a and T b substituted from Eqs (2.38) and (2.39). The same 
conclusion can be drawn with respect to the phase angle shown in fig. 10. 

The spacial variation of the amplitude ratio |T(r)| derived from Eqs. (2.17) and 
(2.40) is plotted in fig. 11. These computations were made a, an angular frequency 
of co/m, =0.1 for the type B thermocouple with the dimensions listed in table 2. It 
is interesting to note the small asymmetry in the amplitude ratio on the left and 
right side of the junction. This asymmetry is the result of the small differences in 
the physical properties across the junction as listed in table 2. 

The steady state amplitude ratio |T(0)| for a type T thermocouple is plotted in 
fig. 12 using equations (2.17) and (2.40). The dimensions and material properties 
are listed m table 2. The average material properties listed in Table 2 were used to 
compute the amplitude ratio of the frequency response with the one material 
solution, Eq. (2.17). Also plotted in fig. ,2 is the amplitude ratio derived from the 

two material solution, Eq. (2.40). In the latter case, the individual material 
properties also listed in table 2 were used. 

As indicated in fig. 12, the amplitude ratio for the approximate one material 
steady-state frequency response of a type T thermocouple is distinctly different 
from the two material approximate solution. This is a consequence of unequal 
material properties across the junction for the type T thermocouple. Also shown 
in fig. 12 is a numerical solution of the second order ordinary differential equation 
for the temperature, Eq. (2.6). The numerical finite difference solution of the 
ndary value problem of Eq. (2.6) matches both the temperature and heat flux at 
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x = 0 and x=±l in the schematic of fig. 2. It is clear from fig. 12 that the 
approximate two material analytical solution, equation (2.40), is in good agreement 
with the numerical results despite some differences at low frequencies 0 / 0 , < .04. 

Similar results are indicated for the phase angle <j>in fig. 13. 

The special variation of the amplitude ratio |T(*)1 derived from Eqs. (2.17) and 
(2.40) is graphed in fig. 14. These computations were made at an angular frequency 
of fti/o), =0.1 for the type T thermocouple. As evident in fig. 14, the two material 
solution derived from Eq. (2.40) accurately represents the features of the 
asymmetry associated with a type T thermocouple. In particular, the relatively 
large resistance to axial heat conduction in the constantan wire on the right of the 
junction is reflected in the larger values of the amplitude ratio PX**. Also shown 
in fig. 14 is the approximate one material solution represented by Eq. (2.17). In the 
latter case, the average values for the materia, properties of a type T thermocouple 
were used as listed in Table 2. Therefore, one can conclude that for thermocouples 
whose material properties on either side of the junction are markedly different, 
the two material solution developed in this paper is a substantial improvement in 
accuracy both in frequency response and in temperature distribution along the 


wire. 


2.5 EXPERIMENTAL PROCEDURE 

In the experiment described below, thermocouple sensors were exposed to a 
constant velocity air stream of varying temperature. In particular, the dynam.c 
response of the thermocouple sensors to a square wave temperature profile was 
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measured for a range of frequencies. 

2.5.1 Rotating Wheel Experiment 

A rotating wheel configuration was used to deliver the test airstream to the 
proposed sensors. A similar experimental apparatus is described in detail by 
Elmore et al. (1986). A schematic of the rotating wheel apparatus used in the 
present experiment is shown in figure 15. As the wheel rotates, the holes pass the 
two air supply tubes which allow slugs of hot and cold air to alternately enter a 
transition tube assembly mounted directly above the rotating wheel. In the 
transition tube the slugs of hot and cold air coalesce into a single air stream 
providing roughly a square wave temperature profile covering a range of 
frequencies from roughly 1 to 30 Hz. 

The analog temperature signal was digitized with a Data Translation DT-2801 
board mounted in an expansion slot of an IBM AT compatible computer as shown 
in figure 15. The ASYST software loaded on the hard disc of the personal 
computer provided a flexible system for data storage, manipulation and display. 

The true temperature profile of the airstream was measured with a constant 
current anemometer (TSI 1054- A) and sensor (1226 PI 2.5). 

2.5.2 Thermocouple Construction 

Thermocouple wire of the desired length and type is threaded into the four 
hole ceramic with the thermocouple end last. Three or four kinks are made in 
each wire near the thermocouple end so the wire must be firmly pulled into the 
tube leaving enough wire sticking out to make the thermocouple. Drops of epoxy 
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are picked up with a piece of .010" diameter wire, added to the ceramic tube at A 
(see fig. 16) and pushed down around both wires. The kinks and epoxy firmly 
fasten the wires in the ceramic so they do not twist when the free ends are 
manipulated for electrical connections. 

The junction at B in figure 16 is made by cutting the large wire about half way 
through with a razor blade, laying the small diameter wire in the cut, and welding 
the cut closed with a stored energy spot welder. After the weld, the excess wire 
ends that protrude through the weld area are bent and broken off to clean up the 
junction. 

For the chromel alumel couples, the junction at C of figure 16 was made using 
a stored energy spot welder. The wires are crossed and welded and then trimmed 
with a razor blade and bent with tweezers under a microscope until collinear. 

A stored energy spot welder will not work for the copper constantan couples. 
For these couples the junctions were made by silver soldering. The silver solder 
wire is coated with flux and the end heated with a torch until a drop forms. A 
twisted pair of thermocouple wires is pushed briefly into the flux coated solder 
drop. The solder will wet the pair up to where the twist stops. Again using a razor 
blade and tweezers under a microscope, the thermocouple is bent and trimmed 
until collinear. 

2.5.3 ASYST Software 

ASYST software was developed to acquire temperature data simultaneously 
from the thermocouple and constant current anemometer. Data were digitized for 
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two channels at a sampling rate of 512 Hz per channel for a total sample time of 
two seconds. The ASYST software code is listed in figure 17. 

The ASYST code first plotted the temperature profiles versus time from the 
thermocouple and anemometer and an example of the plot is shown in figure 18. 
The ASYST software next took the Fast Fourier Transform of the temperature data 
in each channel and recorded the amplitude ratio and phase angle between both 
channels at the first harmonic for the square wave. These data are discussed in the 
next section. 

2.6 EXPERIMENTAL RESULTS 

Initial tests with a signal generator and the output from an RC circuit indicated 
that the ASYST code and data acquisition hardware were operating properly. 
Several test sensors were constructed and tested with the dimensions listed in 
table 3. In each case measurements of the amplitude and phase angle were 
compared with the appropriate theory. A discussion of the results is given below. 

2.6.1 First Order Response 

The amplitude ratio and phase angle were measured with the type K 
thermocouple listed first in table 3. The lengths of the thermocouple and support 
wires in this case were chosen to eliminate the effects of axial heat conduction. 
The experimental data representing the amplitude ratio are plotted in figure 19 
along with the theory representing a first order frequency response. The 
experimental data were correlated with a natural frequency tOj = 5.5 sec -1 and a 
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thermal diffusivity of a = 059 cm 2 / sec. 

As indicated in figure 19 the data has been correlated with the theory to within 
10%. Although the experimental data is reproducible, there is a noticeable drift of 
the data relative to the theory with increasing wheel frequency The drift is 
probably due to a small increase in the natural frequency <w, of the thermocouple 
wire with increasing wheel frequency <y. 

The natural frequency co, of Equation (2.10) depends on accurate predictions of 
the convective heat transfer coefficient h from Equation (2.41). However, there 
may be a small change in either the mean air velocity or unsteadiness in the mean 
velocity which would lead to small changes in the heat transfer coefficient h and 
the natural frequency a ),. 

Also plotted in figure 20 are experimental data representing the phase angle <j> 
for the first order frequency response of the type K thermocouple. In this case the 
data is correlated to within 7% of the theory. The experimental data for the phase 
angle also indicate a small drift in the natural frequency £U t as was apparent in 
figure 19. 

2.6.2 Frequency Response (one material) 

The amplitude ratio and phase angle were measured with the second type K 
thermocouple listed in table 3. The lengths of the thermocouple and support 
wires in this case were chosen (i.e., shortened) to demonstrate the effects of axial 
heat conduction. The experimental data representing the amplitude ratio are 
plotted in figure 21 along with the theoretical prediction given by Eq. (2.32). Both 
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the data and theory were correlated by assuming that the material properties are 
uniform across the thermocouple junction. In this case the data were correlated 

with a natural frequency &>, = 5.0 sec 1 and a thermal diffusivity of 

a = .059 cm 2 / sec. 

As indicated in figure 21 the data has been correlated with the theory to within 
10% over most of the range of wheel frequencies. The data is reproducible but 
again a noticeable drift of the experimental data exists relative to the theory as was 
discussed in the previous section. 

Also plotted in figure 22 are experimental data representing the phase angle <j> 
for the frequency response of the shortened type K thermocouple. The data is 
correlated to within 7% of the theory. The experimental data for the phase angle 
also indicate a small drift in the natural frequency w v 
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2.7 NOMENCLATURE 


A = constant of integration 

B = constant of integration 

c = material specific heat ~ °^ _l ) 

D = thermocouple wire diameter (cm) 

G = \+i(o)/(O n ) 

h = heat transfer coefficient ( J — cm~ 2 - s~ l - “AT -1 ) 

i = unit imaginary number 

k = material thermoconductivity (■/ — cm -1 — s _1 — °K~') 

kf = gas thermoconductivity ( J -cm~ ] - s~ ] - “A" 1 ) 

l = length of small thermocouple wire (cm) 

L = length of large thermocouple wire (cm) 

Nu = Nusselt number (= hD/kf) 

Pr = Prandl number (= v // a ) 

Q = D 2 2 q 2 

M L 

Qt = 

K 103 

Q S = *2<7 4 

q = {Gfyf 

Re = Reynolds number (= vD/ v f ) 

t = time (s) 
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T 


T f 

T g 

To 

— » 

T 

1 to 


V 

X 


steady-state frequency response 
amplitude of periodic gas temperature 
gas temperature (°^) 
mean gas temperature (°^) 

complex amplitude of periodic wire temperature (°^) 

local wire temperature (°K) 

gas velocity (cni-s ') 

axial distance from center of wire (cm) 


Greek Symbols 

a = thermal diffusivity (cm 2 -s ') 

y = a / ft)„ (cm 2 ) 


ft) 

<D 

P 

A, 

A 2 


kinematic viscosity of gas ( cm 2 - j ') 

angular frequency (s ') 

natural frequency of wire ( = 4 /j / pcD){s ') 
phase angle 

-3 

material density (gm - cm ) 
dimensionless function 
dimensionless function 
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2.9 TABLES AND FIGURES 
Table 1 - Properties of One Material Wire (Type B) 

Dimensions (rm) 


a 


Dg 


.025 


.05 


.0076* 


.038 


Average P roperties of Tyg e P 



*co : =32.9 sec -1 


Ug/D, 

2 

5 


I 

.2 

.1 


L 

.35 

.2 


Gas Properties 

T 0 = 900 °K 
M = .26 
A = 1 atm 
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Table 2 - Properties of Two Material Wire 



El 

.0076 


D2 D2ZD1- 
.038 5 


I L 

.1 -2 



P = 1 atm 


V = 50 m / sec 


Pt - 6% Rh 
pt - 30% Rh 
Average 




2.73 


2.86 


2.8 


.238 

.190 

.214 



Copper 

Constantan 

Average 



3.44 116 

3.48 067 

3.46 614 



YJ\re l ocation 
Type B 

pt - 6% Rh c °PP er 

pt - 30% Rh Constantan 
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Type B (DVD. = 5) 



oJ O" 

.is W 
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Amplitude ratio vs. distance for o 
from Eqs. (24) and (31). Crosses are 
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Type B (D 9 /D, = 2) 
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Type B (Dt/D, = 5) 
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Type B (D,/D.= 
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13. Phase angle vs.angular frequency for two material wire. Solid line is Eq. 
(40). Crosses are approximate on material solution Eq. (17). 
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14. Amplitude ratio vs. distance for two material wire. Solid line derived 
from Eqs. (36a) and (36b). Crosses are approximate one material solution 
Eqs. (14) and (16). 


HIGH PRESSURE 
AIR 



15. Rotating wheel experiment. 
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FORGET. ALL 

REAL DIME 1024 ] ARRAY TIME 

COMPLEX DIME 1024 ] ARRAY TRANS 

REAL DIME 6 ] ARRAY ANS 

REAL DIME 1024 ] ARRAY ZMAGO 

REAL DIME 1024 ] ARRAY ZARGO 

INTEGER SCALAR NUM 

REAL SCALAR MEN 

REAL SCALAR CD 

REAL DIME 1024 ] ARRAY FREQS 

REAL DIME 1024 ] ARRAY SIGNAL 

INTEGER dimE 1024 , 2 ] array DATA. BUFFER 


1. CD := 

.0000 DATA. BUFFER := 

LOAD. OVERLAY ACQUIS. SO V 
. 0000 1 . 0000 A/D . TEMPLATE 

DATA . BUFFER TEMPLATE . BUFFER 
CYCLIC 

2048 . 0000 TEMPLATE . REPEAT 
CD CON VERS I ON. DELAY 
DEMO. TEMP LATE A/D. I NIT 
A/D. I N> ARRAY 
LOAD. OVERLAY WAVEOPS.SOV 
1024 REAL RAMP 1 - 2 * CD * 


DEMO. TEMPLATE \ DEFINE AN A/D TEMPLATE 

\ DECLARE ARRAY AS A TEMPLATE BUFFER 
\ SET TEMPLATE BUFFER TO CYCLIC MODE 
\ SET REPETITIONS FOR I/O INSTRUCTION 
\ CONVERSION RATE (MSECS/SAMPLE) 

\ INITIALIZE CURRENT A/D TEMPLATE 
\ A/D INPUT TO TEMPLATES BUFFER 

1000 / TIME := \ SET TIME AXIS 


1 ] 2048 - \ CHANNEL 1 ON STACK (ANEMOMETER) 


\ SMOOTH DATA (CYCLES/POINT) 


DATA. BUFFER XSECTC 
SIGNAL := SIGNAL 

-2048 / 10 * -55.274 * 353. + \ CALIBRATE (DEGREES C) 

SIGNAL := SIGNAL 
0.5 SET . CUTOFF . FREQ 
SMOOTH SIGNAL := SIGNAL 
MEAN MEN := SIGNAL 

MEN - SIGNAL := \ CENTER ON ORIGIN 

TIME SUB[ 1 , 500 , 1 ] SIGNAL SUB[ 1 , 500 , 1 ] 

XY. AUTO. PLOT SIGNAL \ PLOT CHANNEL 1 

FFT \ TAKE FFT 

TRANS := TRANS ZMAG ZMAGO := \ MAGNITUDE OF FFT 

5 SET.#. POINTS 
1 SET # OPTIMA 

ZMAGO SUB[ 1 , 120 , 1 ] LOCAL. MAXIMA \ FIND INDEX AND MAX OF MAGNITUDE 
SWAP NUM ■ = NUM ^ \ INDEX OF MAX MAGNITUDE 

1 - CD / io 24 . / 1000 * PI * \ FREQUENCY AT MAX MAGNITUDE 

ANS [ 1 ] := ANS [ 1 ] DROP ANS [ 2 ] := 

TRANS [ NUM ] DEG ZARG ANS [ 3 ] := \ ARGUMENT AT MAX MAGNITUDE 


DATA. BUFFER XSECTC 2 ) 2048 - 
SIGNAL := SIGNAL 
.124 * SIGNAL := SIGNAL 
0.5 SET. CUTOFF. FREQ 
SMOOTH SIGNAL := SIGNAL 
MEAN MEN := SIGNAL 


\ CHANNEL 2 ON STACK (THERMOCOUPLE) 

\ CALIBRATE (DEGREES F) 

\ SMOOTH DATA (CYCLES/POINT) 


11:42:33.77 


ASYST Version 3.00 

Page 1 TEMP3 . FOR 03/13/91 
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MEN - SIGNAL := 

TIME SUBC 1 , 500 , 1 ] 

SIGNAL SUBC 1 , 500 , 1 ] 

XY. DATA. PLOT SIGNAL 
FFT 

TRANS := TRANS ZMAG ZMAGO := 
5 SET.#. POINTS 
1 SET.#. OPTIMA 


1 - CD / 1024. / 1000 * PI * 

ANS C 4 ] ;= ANS C 4 ] DROP t\ no i o 
TRANS [ NUM ] DEG ZARG ANS [ 6 ] := 

CR 

FREQ MAG PHI " CR 

ANS Cl] ? CR 
ANS C 5 ] ANS C 2 ] / ? CR 
ANS C 3 ] -1 * ANS C 6 ] + 180 PI / 
FORGET. ALL 


\ CENTER ON ORIGIN 


\ PLOT CHANNEL 2 
\ TAKE FFT 
\ MAGNITUDE OF FFT 


AND MAX OF MAGNITUDE 
MAGNITUDE 
MAX MAGNITUDE 

\ ARGUMENT AT MAX MAGNITUDE 


CD * 1000 / ANS C 1 ] * - ? CR 


ZMAGO SUBC 1 , 120 
SWAP NUM := NUM 


] LOCAL. MAXIMA \ FIND INDEX 
\ INDEX OF MAX 
\ FREQUENCY AT 
] • = 


ASYST Version 3.00 

Page 2 TEMP3 . FOR 03/13/91 11:42:34.65 


17. ASYST software 


82 







85 


20. Phase angle vs. angular frequency, tw, = 5.5 sec \ a- .059 cm 2 / sec. 
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3.0 MEASUREMENT OF FREQUENCY RESPONSE 


Experimental measurements are made for the steady-state frequency response 
of a supported thermocouple wire. In particular, the effects of axial heat 
conduction are demonstrated for both a supported one material wire (type K) and a 
two material wire (type T) with unequal material properties across the junction. 
The data for the amplitude ratio and phase angle are correlated to within 10% with 
the theoretical predictions of Forney and Fralick (1991). This is accomplished by 
choosing a natural frequency co n for the wire data to correlate the first order 
response at large gas temperature frequencies. It is found that a large bead size, 
however, will increase the amplitude ratio at low frequencies but decrease the 
natural frequency of the wire. The phase angle data are also distorted for imperfect 
junctions. 


3.1 INTRODUCTION 

The evaluation of jet engine performance and fundamental studies of 
combustion phenomena depend on the measurement of turbulent fluctuating 
temperatures of the gas within the engine (Dils and Follansbee, 1976). Historically, 
these temperatures have been measured with thermocouples. The advantages of 
thermocouples are their low cost, reliability and simplicity since they do not 
require optical access or elaborate support electronics. However, the design of a 
thermocouple represents a compromise between accuracy, ruggedness and rapidity 
of response. 
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For example, the measurement of fluctuating temperatures in the high speed 
exhaust of a gas turbine engine combustor is required to characterized the local gas 
density gradients or convective heat transfer (Fralick, 1985). Although 
thermocouples are suitable for the measurement of high frequency temperature 
fluctuations (<1 KHz) in a flowing gas or liquid, the measured signal must be 
compensated since the frequency of the time dependent fluid temperature is 
normally much higher than the corner frequency of the thermocouple probe 
(Scadron and Warshawsky, 1952). Moreover, the amplitude and phase angle of the 
thermocouple response may be attenuated by axial heat conduction for rugged 
thermocouples of finite length (Elmore, et al.; 1983, 1986). 

In the present study, measurements are made of the steady-state frequency 
response of supported thermocouple wires. In particular, the effects of axial heat 
conduction are demonstrated for both a supported one material wire and a two 
material wire with unequal material properties across the junction. The 
measurements are presented in the form of the amplitude ratio and phase angle 
and these data are compared with the recent analytical results of Forney and 
Fralick, (1991). 

3.2 THERMOCOUPLE FABRICATION 

A supported thermocouple wire shown in figure 1 consists of a small diameter 
thermocouple element (horizontal wire) mounted on two large diameter 
supported wires (vertical wires). The large diameter support wires are mounted in 
a ceramic stick that assumes the mean gas temperature T 0 . Since a supported 
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thermocouple is normally constructed manually, there are a number of physical 
imperfections that can affect the response. For example, three joints must be 
fabricated consisting of the thermocouple junction and joints at both ends of the 
thermocouple element. These three joints (welded or soldered) and the epoxy 
used to secure the support legs in the ceramic stick are darkened in the schematic 
at the bottom of figure 1. As indicated in figure 1, the supported thermocouple 
may not be symmetrical and there may be additional material or beads at the joint 
locations. 

3.2.1. Physical Dimensions 

The thermocouple support wires and the smaller thermocouple elements were 

analyzed by Elmore et al. (1983) for stresses at the maximum flow conditions for 

the FI 00 gas turbine engines of T = 1680 °K, P = 2.01 x 10 6 Newton/meter 2 (202 psia) 

and a gas Mach number M = 0.36. The sensor location for these calculations 

corresponded to a wall position between the vanes downstream of the combustor 

exit. Elmore et al. (1983) also assumed an allowable stress in the thermocouple 
7 2 

wire of 4.39 x 10 Newton/meter (6360 psi). The maximum allowable length-to- 
diameter ratio L/D for the support wire was thus calculated to be 6.5 while the 
maximum value for the thermocouple element was found to be 15.5. 

Elmore et al. (1986) also limited the smallest diameter of the thermocouple 
wire to 76 fim (. 003 in) because of limitations to junction fabrication. These 
constraints limit the length of the thermocouple element on each side of the 
junction to roughly 1 mm or to a total length of 2 mm for the thermocouple wire 


91 



suspended between the support legs. Moreover, if the diameter of the supports is 
380 pm {. 015 in), the length of the support legs are limited to roughly 2 mm. 

3.2.2. Materials 

In the present study thermocouples are constructed of either type K (chromel- 
alumel) or type T (copper-constantan) materials. The properties of these materials 
including the product of specific heat and wire density pc and the thermal 
diffusivity a = k/pc are listed in table 1. The type K thermocouple was chosen 
since one normally assumes that the material properties are equal across the 
junction. In contrast, the type T thermocouple has a distinctly different 
thermoconductivity across the junction. As a result, the temperature profile is 
asymmetric for the type T thermocouple and theoretical expressions derived to 
predict the frequency response must reflect these features. 

The theoretical analysis of a supported thermocouple assumes that the radial 
temperature gradients in a wire cross section are negligable. In this case the heat 
conduction is one - dimensional down the axis of the wire. Schematics of both a 
one and two material supported wire are represented in figure 2. Here, wire 
diameters and lengths are illustrated and their values are listed in table 1. 

The positive lead for both the type K and T thermocouples is positioned on the 
left of the junction in fig. 2 as indicated at the bottom of table 1 . If the material 
properties are nearly equal across the junction, as is the case for a type K 
thermocouple, one can define two natural frequencies for wire regions 1 and 2 in 
fig. 2 where CD, = 4h i jpcD, and co 2 =4h 2 /pcD 2 . For the type T thermocouple, 
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however, there are four natural frequencies. For example, two of these frequencies 
are co, = Ah,/p,c,D x and co 2 = 4h,/p 2 c 2 D, on both sides of the junction. 


3.2.3 Construction 

Thermocouple wire of the desired length and type was threaded into the 
ceramic with the thermocouple end last. Three or four kinks were made in each 
wire near the thermocouple end so the wire must be firmly pulled into the tube 
leaving enough wire sticking ou, to make the thermocouple. Drops of epoxy were 
picked up with a piece of .010 in. diameter wire, added to the ceramic tube at A (see 
fig. 3) and pushed down around both wires. The kinks and epoxy firmly fasten the 


wires in the ceramic so they do 
electrical connections. 


not twist when the free ends are manipulated for 


The junction at B in fig. 3 was made by cutting the large wire about half way 
through with a razor blade, laying the small diameter wire in the cut, and silver 


soldering the cut closed. After the solder, the excess wire ends that protrude 
through the soldered area were bent and broken off to clean up the junction. 

The type T thermocouple junctions a. C in fig. 3 were constructed successfully 
by one of two procedures. In the firs, procedure, the silver solder wire was coated 
with flux and the end heated with a torch until a drop formed. A twisted pair of 
thermocouple wires was pushed briefly into the flux coated solder drop. The 
solder wetted the pair up to where the twist stops. Using a razor blade and tweezers 
under a microscope, the thermocouple was bent and trimmed until collinear. The 
silver solder technique was also used to fabricate the chromel alumel (type K) 
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thermocouples. 

in the second procedure, the thermocouple junctions were fabricated by the 
Research Instrumentation Branch at the NASA Lewis Research Center. In this 
case, the wires were cut with a razor biade to produce a flat edge perpendicular to 
the wire axis. The wires segments were then mounted in a fixture such that the 
two faces of the junction were he.d together by springs. The junction was then 
formed by laser heating. The laser used was a KORAD modei KWD Nd: VAC 
laser, operating at a wave length of 1060 namometers. The power settings depend 
on the wire size and material, bu, for the thermocouples described in this report, a 
pulse duration of 4 ms was used delivering a total energy of approxtmately 2 
joules. One example of a type T junction (D = 76 pm) formed by laser heating are 

shown in fig. 4. 

For the chrome, a.umel couples (type K), the junction a. C in fig. 3 could also be 
fabricated using a stored energy spot welder. In this case, the wires were crossed 
and welded and then trimmed with a razor blade and ben. with tweezers under a 
microscope until collinear. This method proved to be unsuccessful, however, 
since an apparent contact resistance was created by the procedure and the dynamic 
response of the thermocouple was inadequate. 


3.3 EXPERIMENTAL APPARATUS 

in the present experiment, thermocouple sensors are exposed to a constant 
velocity air stream «19 m/sec) of varying temperature. In particular, the dynamic 
response of the thermocouple is measured for a square wave temperature profile 
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of varying frequency. A rotating wheel configuration is used to deliver the test air 
stream to the proposed sensors. A similar experimental apparatus was described 
in detail by Elmore et al. (1986). 

A schematic of the rotating wheel apparatus used in the present experiment is 
shown in fig. 5. As the wheel rotates, holes pass the two air supply tubes (3/4 in. 

ID copper) that allow slugs of hot (^55 °C) and cold (~30°C) air to alternately enter 
a transition tube assembly mounted directly above the rotating wheel. In the 
transition tube the slugs of hot and cold air coalesce into a single air stream 
providing roughly a square wave temperature profile covering a range of gas 
temperature frequencies from roughly 1 to 30 Hz. A 3/4 in. to 1/2 in. smooth 
copper adaptor was inserted at the end of the coalescing tube (location of 
thermocouple) in fig. 5 to increase the air velocity to 18 m/sec. 

The analog temperature signal is digitized with a Data Translation DT-2801 
A/D board mounted in an expansion slot of an IBM AT compatible computer as 
shown in fig. 5. ASYST software was loaded onto the hard disc of the personal 
computer and this provided a flexible system for data storage, manipulation and 
display. The true temperature profile of the airstream is measured with a constant 
current anemometer (TSI 1054- A) and sensor (1226 PI2.5). 

3.4 EXPERIMENTAL PROCEDURE 

3.4.1 Calibration 

The air velocity at the exit of the coalescing tube in fig. 5 (location of 
thermocouple) was maintained at roughly 17.5 to 18 m/sec. The voltage regulator 


was adjusted such that the temperature of the hot line was maintained at an 
elevated temperature of <55 °c. One hole of the rotating wheel was positioned in 
a stationary position over the unheated line. The voltage output from both the 
thermocouple and constant current anemometer were recorded at the input 
terminals to the A/D board. The true temperature was also recorded with a 
thermometer suspended in the high velocity air stream. This procedure was 
repeated for several temperatures with the hole positioned over the hot line 
covering an air temperature range from 30°C to 55°C. It was found that both type 
K and T thermocouples had a constant calibration factor (°C/volt) over the small 
temperature range. The constant current anemometer, however, provided a non- 
linear response with the calibration factor varying with air temperature. Thus, a 
least square curve fit was determined to provide the gas temperature with the 
anemometer output. The calibration factors described above were then installed 
on the ASYST data acquisition software. 

The thermocouple and constant current anemometer are mounted in the 
constant velocity air stream with both sensing elements parallel and separated by 
approximately 1 mm. The period or frequency of the square wave temperature 
profile is controlled over the range from 1.5 to 30 Hz with adjustment of the motor 
connected to the rotating wheel shaft. 

3.4.2 Data Acquisition 

ASYST software was developed to acquire temperature data sequentially from 
the thermocouple and constant current anemometer. The data are digitized for 
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two channels at a sampling rate of 512 Hz per channel for a total sample time of 
two seconds. An example of the gas temperature profile for a square wave 
frequency of 2.5 Hz is illustrated in fig. 6 where the smooth profile is the 
thermocouple response. After data acquisition, the Fast Fourier Transform (FFT) 
is taken for each channel. The largest peak (first harmonic) in the amplitude of 
the FFT is located for each channel which provides the frequency of the 
temperature profile. The ASYST software also records the amplitude ratio and 
phase angle between both channels at the first harmonic. In the latter case, the 
software corrects the phase angle for the conversion delay between channels. The 
amplitude ratio and phase angle are recorded for each frequency setting of the 
rotating wheel. 

3.4.3 Data Correlation 

Once the amplitude ratio and phase angle are recorded for a large number of 
gas temperature frequencies co for a given thermocouple, the amplitude ratio and 
phase angle are plotted on separate graphs. Both quantities are plotted as a 
function of the ratio of the angular frequency of the square wave temperature 
profile co to the natural frequency a>, of the positive thermocouple element. The 
natural frequency co^ = 4h l /p ] C i D ] is chosen to fit the amplitude ratio at large 
frequencies since the amplitude ratio should coincide with the first order response 
(independent of axial heat conduction) at large q). 

The lengths of the thermocouple element 1 and support wire L are then chosen 
such that theoretical predictions of the amplitude ratio are correlated with the 
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data at low gas temperature frequency co- The remaining parameters representing 
wire diameters and material properties were recorded from the literature and the 
manufacture s specifications. These values were also introduced into the analytical 
expressions used to predict both the amplitude ratio and phase angle. 

Theoretical predictions of the amplitude ratio and phase angle are plotted along 
with the data. For the type K thermocouples the material properties are assumed 
to be equal across the junction with the average values of the thermal diffusivity 
a and the product of the density and heat capacity pc listed in table 1. Predictions 
of the amplitude and phase angle for the type K thermocouple were derived from 
the expression representing the frequency response (Forney and Fralick, 1991) 


7(0) 


-^-(cosht/ 2 L- l)--t-cosh<? 2 L 

^2 fdj + _L 

cosh q x l cosh q 2 L + ^sinh q x l sinh( 7 2 L G, 


(3.1) 


where 


Q = 


■D 2 *7 2 


and 


^ = 7T’ G a =\ + i((o/o) n ), co n = 

V Yn ton 


4/*„ 

pcD n 


Here, the subscripts 1 and 2 refer to the wire elements 
schematic of figure 2 and T(0) = |T(0)| exp(t’O) where |T(0)| is 
and is the phase angle. 

Theoretical predictions of the amplitude ratio and phase 
thermocouple were derived from the expression 


in the one-material 
the amplitude ratio 

angle for the type T 
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7(0) = — < (9, sinh < 7 2 / T a + — (cosh#,/- 1 ) + sinh < 7 ,/ T b +— (cosh<? 2 / - 1) 

A i [ G i J I G 2 


where 


0 _ K2i 
k 2 q 2 


= Q , cosh qj sinh q 2 l + cosh^/sinhg,/ . 

The boundary temperatures T a , Tb in Eq. (3.2) are given by the expressions 


T = — + jsinh^L \ ~ -J- ]cosh^ 4 (/ + L)- -J- 
G, \ g a g J G 4 


— — [^ cosh ^ 3 / sinh g 4 (/ + L) + coshc/ 4 (/ + L) sinh ^ 3 /] if — 


T. = — + |< 2 f sinh< 7 4 L jcosh <y 3 (/ + O --7 

b g a G 3 GJ G 3 


1 1 

[cosh^ 4 /sinh^ 3 (/ + L) + 2 J cosh^ 3 (/ + L)sinh< 7 4 /] f — 

Ga V A, 


where 
O = kl<h 

k 2 q, 

A 2 = Q s cosh^ 3 (/ + L)sinh^ 4 (/ + L) + cosh^ 4 (/ + L)sinh q 3 (l + L) 

Here, the wire elements 1 to 4 are shown in the schematic of a two-material 


thermocouple wire shown in figure 2 . 


99 



3.5 RESULTS AND DISCUSSION 

A type K thermocouple was silver soldered with a thermocouple element 
length of 1 = 0.9mm and a support wire length of L = 2.2mm. The thermocouple 
and support wires had diameters of D, = 76 f±m and D 2 = 380/im, respectively. The 
amplitude ratio and phase angle were measured over a range of temperature 
square wave frequencies of 1.3 to 35 Hz and these data are plotted in figure 7. The 
data are also compared with the theoretical predictions derived from Equation (1). 

A value of the natural frequency ft), was chosen to correlate the data when 
co/ft), > 1. The value of ft), chosen was within 10% of the predicted value based on 
the air velocity and thermocouple wire diameter. As indicated in figure 7, all of 
the data are correlated with the theory to within 10%. It should also be noted that 
a small bead did exist at the junction with a diameter less than twice the diameter 
of the thermocouple element wire. Moreover, the attenuation of the signal 
because of axial heat conduction is properly accounted for with the theory 
provided the wire length 1 is increased in Eq. (3.1) by 20% to account for the bead at 
the junction. 

A type T thermocouple junction was purchased from Omega Engineering and 
the thermocouple wire was silver soldered at the support legs. This thermocouple 
had an element length of 1 = 1.25mm and a support wire length of 2.2mm. The 
thermocouple and support wires had diameters of £>, = 16fim and D 2 = 380/t/n, 
respectively. The amplitude ratio and phase angle were measured over a range of 
temperature square wave frequencies of 1.3 to 30 Hz and these data are plotted in 
figure 8. The data are also compared with the theoretical predictions derived from 
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Equation (3.2). 


A value of the natural frequency to, =20 sec" „ as chosen to correlate the data 
for the amplitude ratio at large frequencies oj/oj, > 1. As indicated in figure 8, the 
data for the amplitude ratio drift below the theoretical prediction by 10% for 
increasing gas temperature angular frequencies «. Although the effect of axial 
heat conduction is properly accounted for in figure 8, the phase angle data also 
appear to be distorted and to he significantly below the predictions. These effects 
are apparently due to the large bead diameter at the junction (greater than twice 
the diameter of the thermocouple wire). Here, as before, it was necessary to 
increase the wire length in Eq. (3.2) by 20% to correlate the data. 

A type T thermocouple junction without a bead was constructed by laser 
heating with the support legs silver soldered, as before. The thermocouple had an 
element length of 1-1.25 mm and a support wire length of 2.5 mm. The 
thermocouple and support wires had diameters of D, = 76 pm and D 2 = 380 pm, 

respectively. The amplitude ratio and phase angle were measured and the results 
were compared with theory. 

A value of the natural frequency to, =28 sec' was necessary to correlate the 
amplitude data at large frequencies. This represents a larger air velocity compared 
to the data in figure 8. As shown in figure 9, there appears to be less drift in the 
amphtude ratio below the theory compared to the previous results taken with a 
beaded junction. An improved correlation is also observed for the phase angle 
data plotted in figure 9. The predicted phase angle from Equation (3.2) correlates 
the data well and there appears to be less distortion of the phase angle data 
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compared to the previous results with a bead at the juction. Moreover, no 
adjustments were necessary in the theory to correlate the data. 

Finally, a type T thermocouple junction was constructed with silver solder on a 
twisted pair of thermocouple wires. The junction consisted of a twisted pair of 
soldered thermocouple wires 1.0 mm in length that was bent perpendicular to the 
axis of the wire elements suspended between the support legs. The wire element 
length was 1 - 1.3 mm and the length of the support leg was L = 2.2 mm, as before. 
This thermocouple configurat.on was constructed to determine the effect of an 

extremely large beas size on the frequency response. 

The amplitude ratio and phase angle for the beaded type T thermocouple are 
plotted in figure 10. As indicated in figure 10, the amplitude is larger than 
predicted at low frequencies since the extended exterior surface of the bead (fin) 
enhances the rate of heat transfer to the junction. At large frequencies, >1 

however, the smaller area - to - volume ratio decreases the natural frequency of 
the bead configuration relative to a beadless cylinder and the amplitude data fall 
below the theory. More difficult to describe is the distortion in the phase angle 
data. As indicated in figure 10, the phase angle is shifted to larger values relative 

to the theory. 

3.6 CONCLUSIONS 

Supported type K and type T thermocouple wires were constructed. The 

thermocouple element (D =76^) and wire support (D = 380,1m) lengths were 
approximately 1.0 and 2.0 mm, respectively. Measurements of the amplitude ratio 
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and phase angle that reflect the effects of axial heat conduction are correlated to 
within 10% with the theoretical predictions of Forney and Fralick, (1991). To 
correlate the experimental data, it is necessary to choose a natural frequency co n 
(depends on gas properties and velocity) for the amplitude data at large gas 
temperature frequencies. If there is a bead at the junction (typically twice the wire 
diameter), it is necesary to increase the wire length 1 by 20% in the theory to 
correlate the data. 

The use of silver solder for the junction and wire supports is adequate for both 
the type K and type T thermocouples but will leave a small bead. A spot welding 
procedure for the type K thermocouple, however, appears to introduce a large 
contact resistance since the measured amplitude ratio in this case was only a 
fraction of its predicted value. The use of laser heating to form the junction, 
however, left no bead and required no adjustments to the wire length in the 
theory to correlate the data. Moreover, the use of laser heating for the purpose of 
fabricating thermocouple junctions without beads appears to improve the 
correlation of both the amplitude ratio and phase angle data with theory. In 
particular, the phase angle data from a beadless junction exhibited far less 
distortion compared to theoretical predictions. Also, the wire lengths used in the 
theory to correlate the data were realistic for the laser heated junctions. 

It is found that a large bead at the thermocouple junction increases the 
amplitude ratio at low frequencies but decreases the natural frequency of the wire. 
The phase angle data are also distorted for imperfect junctions. 
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3.7 NOMENCLATURE 


c 

D 

G 

h 

i 

k 

1 

L 


Q 


Qt 


Q s 


q 

T 

x 


material specific heat {J ~ 8 m ' ) 

thermocouple wire diameter (cm) 

l + i(co/co n ) 


heat transfer coefficient (J -cm 2 -s l -°K ') 

unit imaginary number ( = V-l) 

material thermoconductivity ( J-cm ’-5 l -°K ') 
length of small thermocouple element wire (cm) 
length of large thermocouple support wire (cm) 


Ek l 

M L 
k 2*?2 

KQz 

k 2 q, 

[Girt 

steady-state frequency response 

axial distance from junction of wire (cm) 


Greek Symbols 

a = thermal diffusivity (cm 2 -s') 

y = a/(0 n (cm 2 ) 

c 0 = angular frequency ($ ') 
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«>n 

O 

P 

A, 

^2 


natural frequency of wire (= 4h / pcD)(s ') 
phase angle 

material density (gm/cm 3 ) 
dimensionless function 
dimensionless function 
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3.9 TABLES AND FIGURES 


Table 1 Properties of Thermocouples 

Dimensions (cm) 


a 

02 

O 2 /D 1 

1 

.0076 

.038 

5 

.1-.15 


L 


.17-. 22 


Chromel 

Alumel 

Average 


Properties of Type K 


pc 


J 


cm 2 


a 


cm 3 - K _ 


sec 


3.80 

4.10 

3.95 


.048 

.070 

.059 


Properties of Type T 


cm 


pc 

cm 3 - K _ 

a 

sec 

Copper 

3.44 

1.16 

Constantan 

3.48 

.067 

Average 

3.46 

.614 



Wire Location 


Region 

Type K 

Type T 

left 

Chromel 

Copper 

right 

Alumel 

Constantan 
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SUPPORTED THERMOCOUPLE 


IDEAL 


/ 



JUNCTION 



1. Schematic of Supported Thermocouple 
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Fluid Flow 



Fluid Flow 



2. Schematic of One or Two Material Supported Wires 


O' 




4. Photograph of a type T junction formed by laser heating. 
76pm. 


Diameter of wire is 
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HIGH PRESSURE 
AIR 



5. Schematic of Rotating Wheel Apparatus 
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